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ABSTRACT 



We report the detection of the Doppler shadow of the transiting hot Jupiter CoRoT-1 lb. Our analysis is based on line-profile tomog- 
raphy of time-series, Keck/HIRES high-resolution spectra acquired during the transit of the planet. We measured a sky-projected, 
spin-orbit angle A = 0.1° ± 2.6°, which is consistent with a very low-inclined orbit with respect to the stellar rotation axis. We re- 
fined the physical parameters of the system using a Markov chain Monte Carlo simultaneous fitting of the available photometric and 
spectroscopic data. An analysis of the tidal evolution of the system shows how the currently measured obliquity and its uncertainty 
translate into an initial absolute value of less than about 10° on the zero-age main sequence, for an expected average modified tidal 
quality factor of the star (<2'*)~4x 10 6 . This is indicative of an inward migration scenario that would not have perturbed the pri- 
mordial low obliquity of CoRoT-llb. Taking into account the effective temperature and mass of the planet host star (T cff = 6440K, 
Mi, = 1.23 Ma), the system can be considered a new telling exception to the recently proposed trend, according to which relatively hot 
and massive stars (T eff > 6250 K, M+ > 1.2 M Q ) seem to be preferentially orbited by hot Jupiters with high obliquity. 

Key words, planetary systems - planet-star interaction - stars: individual: CoRoT-1 1 



1. Introduction 

Time-series, high-resolution spectroscopic observations of plan- 
etary t ransits allow us to detect the Ro ssiter-McLaughlin (RM) 
effect dRossiterfl92 4: McLaughlin 1924ft and to measure the sky- 
projected system obliquity A, i.e., the angle in the plane of the 
sky between the projections of the planet's orbital angular mo- 
mentum and the star's rotation spin. When a planet transits in 
front of its host star, it subsequently occults different parts of the 
stellar disc. For a rotating star, this results in a distortion of the 
stellar line profiles that change during the transit. If stellar rota- 
tion is the dominant source of line-broadening, high-resolution 
spectroscopy can resolve the stellar line profiles and shows the 
distortion (bump) moving across the line profiles, which is the 
result of the missing starlight bloc k ed by the tra nsiting planet 
(e.g. JCollier Cameron et alBOlOallbtlMiller et al.ll2010l) . On the 
other hand, if stellar rotation is not the dominant broadening 
mechanism, the RM effect is de tected as an anomalous Doppler 
shift of the stellar lines (e.g.. | Queloz et all l2000t Ohtaetall 



l2005HGimenezll2006HGaudi & Wlndl2007tlWinnll2010h~ 

Giant planets with a semi-major axis a < 0.1 AU (hot 
Jupiters) are thought to have formed at large orbital distances 
from their parent star and then migrated inwards to their cur- 
rent position through a) tidal interaction with the proto-planetary 



disc in the pre-main seque nce phase of the star/planet evo- 
lution (e.g., Lin et al.l fl996); b) planet-planet scattering (e.g., 
iRasio&Fordl 119961) : c) gravitational interaction with a third 
outer planet or stellar companion (Kozai mechanism; e.g. [K ozai 
| 1962[ iFabrvckv & Tremainell2007t IWuet al.l 120071: iNaoz et al. 
1201 ll) . The migration mechanism is not completely understood 
yet and the proposed theories predict different final obliquities 
with low values in the case of a migration governed by planet- 
disc interaction and significant misalignment in the case of th e 
other mechanisms ( Triau d et al.l2010t iMorton & Johnsonl201 ll) . 
Measuring the spin-orbit obliquity of hot Jupiters can thus pro- 
vide insights into their migration mechanism and help to discern 
bet ween the riv aling theories. 

IWinn et ail d2010h empirically found that hot Jupiters orbit- 
ing stars with an effective temperature r e ff > 6250 K and a mass 
M+ > 1 . 2Mp) tend to be mi saligned. A similar result was ob- 
tained b v I S chlauf man! d20 10). who found evidence of spin-orbit 
misalignment along the line of sight in system s with massive 
stars and planets. According to lWinn et al.l (120 10). either a planet 
formation scenario dependent on the stellar mass, or an efficient 
tidal realignment of hot Jupiters around cool stars, mig ht account 
for the suspected trend. More recently, Tria udl(l201 ll) suggested 
that the relatively rapid evolution of stars with M* > 1.2 M@ 
combined with an observational bias, might explain the lack of 
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aligned systems for stars with > 6250 K. However, the num- 
ber of exoplanets with secure measurements of the A angle is 
not yet statisti cally significant enou gh to confirm this trend, as 
pointed out by Mo utou et all (1201 lh . 

iGandolfi et al.l ([2010) announced the discovery of the transit- 
ing hot Jupiter CoRoT-1 lb, a 2.33 Mj up planet in a 2.99 days or- 
bit around a T e g =6440 K dwarf star (SpT=F6 V) with a relatively 
high projected rotational velocity (v sin z* = 40 km s _1 ). With the 
aim of detecting the RM effect of the planet, the authors ac- 
quired radial velocity (RV) measurements during the transit, us- 
ing the High Resolution Echelle Spectrometer (HIRES) mounted 
on the Keck I 10 m telescope, at the Keck Observatory (Mauna 
Kea, Hawai'i). Unfortunately, the observations were scheduled 
according to an old, slightly incorrect transit ephemeris, result- 
ing in a partial coverage of the event. Although the HIR ES RV 
data c learly show the RM anomaly of CoRoT-1 1. IGandolfi et alJ 
d2010i) were unable to place strong constraints on the orbit obliq- 
uity, owing to the partial coverage of the transit event and the 
relatively high vsinz* of the star, which strongly affected the 
precision of their RV measurements. 

In this letter, we show how the problems encountered in de- 
riving accurate RV measurements of rapidly rotating stars to 
assess the shape of their RM anomaly can be overcome with 
Doppler imaging techniques, even in the occurrence of a par- 
tial spectroscopic coverage of the transit event. On the basis of 
line-profile tomography applied to the above-mentioned HIRES 
spectra, we succeed in detecting the Doppler shadow of CoRoT- 
1 lb and measuring the sky-projected obliquity of the system. A 
simultaneous reanalysis of the available photometric and spec- 
troscopic data using a Markov chain Monte Carlo (MCMC) ap- 
proach, enabled us to refine the system's parameters. 



2. Observations and data analysis 

As part of NASA's key science project to support the CoRoT 
mission, high-resolution (R « 50 000) tran sit spectroscopy w as 
carried out using the HIRES spectrograph (Vogt et al.ll 19941) . A 
detailed description of t he adopted inst r umen t set-up and data 
reduction was reported in Gandol fTet al.l(l2010l) . Twelve consec- 
utive spectra of 900 seconds each were gathered during the tran- 
sit that occurred on the night of 1 July 2009. The observations 
started about 1.2 hours before the beginning of the transit and 
erroneously stopped ~45 minutes before the end of the event, 
owing to an incorrect transit ephemeris at the time of the ob- 
servations. The altitude of CoRoT-1 1 increased during the time- 
series observations, with an airmass varying from 1.80 to 1.05. 
Two more out-of-transit spectra were secured on the night be- 
fore and one on the same night, to properly combine the HIRES 
RV measurements with data previously taken with other instru- 
ments. 

The HIRES spectra were acquired using an I2 absorption cell 
to track any instrumental drifts and improve the RV precision. 
Since our analysis is based on Doppler imaging of the stellar line 
profile, we used only the echelle orders that were not contami- 
nated by the I2 absorption lines. The first spectrum was discarded 
as affected by the high airmass of the exposure. This resulted in 
a set of eleven consecutive HIRES spectra covering the wave- 
length range 3800 - 4860 A, with a signal-to-noise ratio ranging 
from 40 to 70 per pixel at 4300 A. 

For each exposure, a composite stellar line profile was com- 
puted using the lea st squares deconvolution (LSD) algorithm of 
iDonati et alJ JT997). This is the line profile that, when convolved 
with a mask of delta functions at the wavelengths of known 
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Fig. 1. Upper panel: Least squares, deconvolved line profiles 
of CoRoT-11 (points with error bars), as obtained from the 
time-series HIRES spectra acquired during the transit of the 
planet. The limb-darkened, rotationally-broadened line profile 
of the best fitting model is shown with solid lines. Radial ve- 
locity (wavelength) increases from left to right, time from bot- 
tom to top. The planet signature is the bump that migrates from 
about -20 kms _I near the middle of the time series, to about 
+ 15 kms 1 in the upper line profile. Lower panel: Residuals 
of the composite spectral line profiles of CoRoT-11 following 
the subtraction of the above-mentioned synthetic profile. Orbital 
phase increases vertically. The two vertical dashed lines mark the 
±v sin i-t, of the star. The three crosses denotes the first, second, 
and third planet's contacts. A moderately bright feature, "travel- 
ling" from the bottom left to the top right corner, appears in the 
residual map as the result of the transit of the planet across the 
stellar disc. 



lines, with relative depths compu ted for an F6 V spectrum us- 
ing a Kumcz atmosphere model (lKurucjfl979h . yields an op- 
timal inverse-variance-weighted fit to the observed spectral or- 
ders. The stellar line mask included 1 140 lines in the 15 HIRES 
echelle orders spanning the wavelength range from 4016 A to 
4856 A. The deconvolved profiles were computed on a linear ve- 
locity scale in the barycentric frame, with a velocity increment 
of 1.5 kms^'pixeL 1 (FigureQ] upper panel). 

The lower panel of Figure Q] shows the time-series residuals 
of the composite spectral line profile in grey-scale form, follow- 
ing subtraction of the the limb-darkened, rotationally-broadened 
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line profile of the best-fitting stellar model. The trailed spectrum 
reveals a narrow, moderately bright feature moving with a con- 
stant radial acceleration through the composite stellar line pro- 
file from the bottom left to the top right corner. This feature is 
identifiable as the spectral signature of the photospheric light 
blocked by CoRoT-llb as it transits in front of the stellar disc. 
Important results can be drawn from a qualitative inspection of 
the plot in Figure Q] First, the planet's Doppler shadow appears 
over the approaching limb of the star and moves towards the re- 
ceding half, confirming the prograde orbit s cenario suggested b y 
the shape of the observed RV anomaly (iGandolfi et alJ uOlO). 
The absolute value of the radial velocity of the light blocked by 
the planet at the beginning of the transit is ~ 20 kms -1 , which 
agrees with the expected value for a planetary transit with an 
impact parameter b = 0.81 and for a stellar projected rotational 
velocity v sin /* = 40 km s _1 (see TableQ]). At the phase of mid- 
transit, the trail has a zero barycentric velocity relative to the 
star, i.e., the planet crosses the sky-projected stellar spin axis at 
the mid-point of the transit. Taking into account the non-central 
planetary transit (b + 0), this indicates that the CoRoT-1 1 system 
is closely aligned in the plane of the sky (A ^ 0°). 

The time-series spectra were modelled together with the 
CoRoT photometry and the out-of-transit RV dat a from HIRES, 
HARP S, SOPHIE, and TLS/Coude, as listed in IGandolfi etalJ 
(2010). The system para meters were computed using th e MCMC 
approach described bv lCollier Cameron et al.l (12010a). in which 
the photometry, the out-of-transit RV data, and the least squares 
deconvolved rotation profiles are fitted simultaneously. 

The photometric t ransit model employs the algorithm of 
Mandel & AgoJ (120021) in the small-planet approximation, using 
a four-coefficient non-linear limb-darke ning model w ith coeffi- 
cients interpolated from the tabulation of Claret (2004) to match 
the CoRoT bandpass for a star with the same effective temper- 
ature, metallicity, and surface gravity as CoRoT-11. The ligh t 
contamination fraction of 13.0 + 1.5 % (Gand olfi et al.l [2010) 
was taken into account in modelling th e CoRoT photometry. 
The m odell ing procedure i s descri bed by ICollier Cameron et al] 
d2007l) and iPollacco etaP d2008l) . and yields five parameters: 
the orbital period P, the epoch 7q of mid-transit, the planet- 
to-star area ratio (R p /R+) 2 , the scaled stellar radius R+/a as a 
fraction of the orbit's semi-major axis a, and the dimensionless 
impact parameter b = a cos i p /R+, where i p is the inclination 
of the planet's orbital axis to the line of sight. The stellar mass 
is estimated at each time step using the empirical calibration of 
lEnoch et"aT1 (|2010[). based o n the eclipsing-binary data compi- 
lation of Torres et al. (2010). The out-of-transit RV curve was 
fitted using a circular orbit model of velocity amplitude K. A 
separate zero-point velocity was fitted for each of the different 
RV data sets. 

For the line-profile analysis during the transit, we used 
the line -profile decomposition model of ICollier Cameron et al.l 
(1201 Pal) , in which the stellar profile is modelled as the convolu- 
tion of a narrow Gaussian of width vfwhm (representing the com- 
bined instrumental and intrinsic stellar photospheric line profile) 
and a limb-darkened rotation profile of width v sin /+. The light 
blocked by the planet was treated as a Gaussian of the same 
width as the intri nsic profile, whose stre ngth was modulated us- 
ing the model of Mandel & Agol (2002). The trajectory of the 
planet's signature in velocity space is governed by the projected 
spin-orbit misalignment angle A, the stellar vsin/*, the stellar 
radial acceleration, the transit impact parameter b, and the zero- 
point of the velocity scale on the LSD profile (vo). The latter 
parameter may differ significantly from the systemic velocity 
y established by the out-of-transit RV measurements, as it de- 
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Fig. 2. Correlation diagrams for the joint posterior probabil- 
ity distributions of the four tomography parameters, namely A, 
v sin vfwhm, and vq. The colours of the dots reflect the A% 2 of 
each point, with thresholds at 2.30, 6.17, and 11.8, to show the 
two-parameter one, two, and three-sigma confidence regions in 
yellow, red, and blue, respectively (see the online edition of the 
journal for a color version of this figure). 



pends on the form of the stellar spectrum and the line weights 
and wavelengths used in the LSD procedure. In all, the trailed 
spectrogram was modelled as a function of four new MCMC pa- 
rameters (A, v sin i+, vfwhm, and Vo), in addition to the RV curve 
semi-amplitude K (which is also constrained by the RV solu- 
tion) and the impact parameter b (which is also determined by 
the photometry). Figure |2]shows the correlation diagrams for the 
joint posterior probability distributions of the four tomography 
parameters, as derived from our MCMC analysis. The derived 
system parameters and their uncertainties are listed in TableQ] 

3. Discussion 



After HD189733b dCoIjler Cameron et a l l |2010ab, WASP-33b 
(ICollier Cameron et al.l l2010bl) . and WASP-3b dMiller et al.l 
2010), CoRoT-llb is the fourth transiting exoplanet for which 
the A angle has been derived using a Doppler imaging tech- 
nique. This letter has demonstrated once more the advantage of 
this method relative to standard RV approaches for the analy- 
sis of the RM anomaly, especially in cases where accurate RV 
measurements are rendered unfeasible by the rapid stellar ro- 
tation rate. The RM eff ect of CoRoT-llb was first observed 
by IGandolfi et al.l {2010), who detected the RV anomaly dur- 
ing the transit of the planet. Although the shape of the RV 
anomaly was clearly indicative of a prograde orbit, the sky- 
projected obliquity of the system could not be co nstrained from 
the RV data. The measurements presented in IGandolfi et al.l 
(120101) did not cover the entire transit and were strongly af- 
fected by the high v sin i+ of the star, preventing any accurate 
modelling of the RV anomaly. We reanalysed the existing on- 
transit sp ectra with the line-profile tom ography technique de- 
scribed in Collier Camer on et al. (2010b) and confirmed a pro- 
grade orbit for CoRoT-llb. By tracking the trajectory of the 
missing starlight across the line profile, we measured the sky- 
projected system's obliquity A = 0.1° + 2.6°. The advantage of 
tracking the Doppler shadow of a transiting planet compared to 
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Fig. 3. vsin/+ distributio n of the 18 CoRo T-1 1-like stars from Fig. 4. Forward tidal evolution of the obliquity iff for 



the compilation of lValent i & Fischei] d2005l) . 



4 X 10 6 (solid lines) and (Q ;> = 2 x 10 7 (dotted lines). The grey 
shaded area indicates the age of the star (2+1 Gyr). 



modelling the RV anomaly occurring during the RM effect, is 
that the planet's trail is a straight line and its position can be 
extrapolated even with a partial spectroscopic coverage of the 
transit. Modelling the RV anomaly can lead to a degeneracy of 
the solution, especially if the RV data are noisy and do not cover 
the entire event. 

A combined MCMC analysis of the available CoRoT pho- 
tometry and spectroscopic data, yields the system parameters 
listed in Table Q] A comparison with the previous values re- 
ported in the disco very paper and those published recently by 
Sout hworthl (1201 ll) . illustrates the very good agreement (within 
lcr) of the three sets of results (TableQ}. 

The RM effect allows us to measure only the sky-projected 
angle A. The true obliquity iff between the stellar spin axis and 
the orbital angular momentum can be derived only if one also 
knows the inclinations of the planetary orbit i p and stellar spin 
axis ii, with respect to the line of sight. While the transit light- 
curve modelling provides the former angle, the inclination i+ can 
be derived from the projected rotational velocity v sin /+ and stel- 
lar radius R*, once the rotation period of the star P rol is known. 
Unfortunately, no evidence of spot-induced rotational modula- 
tion has been found in the CoRoT light cu rve, which is ex- 
pected for a F6V star (iGandoffi et al J 120101) . The low activity 
level of CoRoT-1 1 is further confirmed by the lack of detectable 
spot-induced signatures in the composite stellar line profiles (cf. 
Figure [TJ. 

We searched for possible evidence of an equator-on view 
of CoRoT-11 by comparing its vsinz* with the projected ro- 
tational velocity of a sample of randomly oriented stars sim- 
ila r to CoRoT-11. F o r this purpose, we used the compilation 
of I Valenti & Fischer! (|2005), which provides the projected ro- 
tational velocity for 1040 field stars. Figure [3] shows the vsini* 
distribution of the sub-sample of 18 F-type stars, whose parame- 
ters resemble those of CoRoT-1 1. Although affected by a small- 
number statistics, the histogram seems to be asymmetric and 
peaks between 10 kms" 1 and 15 kms -1 . With a projected ro- 
tational velocity of 38.47 + 0.07 kms -1 , CoRoT-11 lies in the 
very right tail of the distribution, suggesting that i+ is likely to 
be close to 90°. 

A constraint on the true obliquity iff can be obtained from the 
variation in the duration of the transits of CoRoT-11. Since the 
star is rapidly rotating, its sizeable quadrupole moment induces 
a precession of the lines of the nodes of the orbital plane on a 



time scale as short as ~ 7 x 10 4 years (lLanza et al. I l201lh . This 
produces a variation in the duration of the transits that can be 
used to constrain the obliquity. The effect is significant even for 
small obliquities if/ < 10° and i s worth searching for in the near 
future (iDamiani & Lanzall201 ll) . 

Our determination of A is consistent with the system cur- 
rently being nearly aligned in the plane of the sky. This makes 
CoRoT-11 an interesting exception to the empirical rule stating 
that st ars with T eS > 62 5 K and M* ^ L2M B have high obliq- 
uities dWinn et al.ll2010h . lLanza et all (1201 ll) pointed out that the 
tidal interaction time scale of the system is comparable to the 
estimated age of CoRoT-11 (2 + 1 Gyr) if the average modified 
tidal quality factor of the star is 4 x 10 6 < (Q ^ 2 x 10 7 . One 
might thus wonder whether the orbit of the planet was initially 
highly tilted with respect to the stellar spin axis and the strong 
tidal torques then gradually aligned the system, damping down 
any initial misalignment. 

To answer this question, we studied the tidal evolution of the 
obliquity, followin g the initial condition scenario described in 
lLanza et aD (1201 ll) . According to this scenario, on the zero-age 
main sequence (ZAMS) the planet would have had an initial or- 
bital period that was nearly synchronous with the stellar spin. 
We assumed that /* = i p = 83.41° (i.e., if/ « A), which implies 
a current value of the stellar rotation period of P rot = 1.74 days. 
Figure|4]shows the forward ti dal evolution of the true obliquity if/ 
in the (Q *) domain found bv lLanza et al.l (1201 ll) for the CoRoT- 
1 1 system, namely, 4xl0 6 $(g'*)~2xl0 7 . With a system age 
of 2 + 1 Gyr, the current A = 0.1° + 2.6° translates into an initial 
obliquity \if/\ < 10° on the ZAMS. We remind the reader that 
(Q *) ~ 4 X 10 6 would imply a non-synchronous initi al state of 
the sy stem, as well as an age younger than 1 Gyr (lLanza et al.l 
1201 ll) . This is in contrast to the estimated age (2+1 Gyr). On 
the other hand, if the star-planet tidal interaction were weak, i.e., 
(6*) ~ 2 x 10 7 , the system would basically not have evolved 
during its lifetime and the current low obliquity would be close 
to that of its initial state on the ZAMS. The possibility that (2 '*) 
might be as small as 10 6 in the case of obliquity evolution, while 
of the ord er of 10 7 - T O 9 in the case of planet migration was sug- 
gested bv lLail (1201 2l) and leads to results similar to those shown 
in Figure [4] pointing again towards an initially low obliquity. 
Finally, an initially high eccentric orbit of CoRoT-llb would 
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Table 1. CoRoT-1 1 system parameters as derived from our MCMC analysis and comparison with the values reported in the literature. 



Parameter 


Symbol 


Gandolfi et al. (2010) 


Southworth (2011) 


This work 


Unit 


Orbital period 


P 


2.994330 ±0.000011 


- 


2.994325 ± 0.000021 


days 


Transit epoch 


T 


2454597.6790 ± 0.0003 


- 


2454690.49838 ± 0.00016 


days 


Transit duration 




2.5009 ±0.0144 


- 


2.44488 ± 0.01800 


hours 


Impact parameter 


b 


0.8180 ±0.0080 


- 


0.8108 ±0.0077 


R* 


Planet-to-star area ratio 




0.011449 ±0.000107 


0.011381 ±0.000483 


0.011600 ±0.000230 




Scaled star radius 


RJa 


0.1451 ±0.0017 


0.1452 ±0.0022 


0.1416 ±0.0024 




Orbit inclination 


b 


83.17 ±0.15 


83.13 ±0.19 


83.41 ±0.17 


deg 


Orbit semi-major axis 


a 


0.0436 ± 0.005 


0.0440 ± 0.0016 


0.04351 ±0.00036 


AU 


Orbit eccentricity 


e 


(fixed) 


(fixed) 


(fixed) 




RV semi-amplitude 


K 


280.0 ± 40.0 


- 


304.0 ± 32.0 


m s _1 


Systemic radial velocity" 


7 


-1.336 ±0.044 


- 


-1.342 ±0.041 


kms" 1 


Star mass 


M+ 


1.27 ±0.05 


1.26 ±0.14 


1.23 ±0.03 


M 


O l III 1 clUl 


R 


1 37 + 03 


1 T.HA + fl flf, 1 

1 . J> I *■+ \J .\J\J L 


1 33 + 04 




Star mean density 


P* 


0.69 ± 0.02 


0.69 ± 0.03 


0.74 ± 0.04 


gem -3 


Star projected rotational velocity 


v sin i+ 


40 ±5 




38.47 ± 0.07 


kms -1 


Sky-projected spin-orbit angle 


A 






0.1 ±2.6 


deg 


Planet mass 


M p 


2.33 ±0.34 


2.34 ±0.39 


2.49 ± 0.27 


Mj up 


Planet radius 


Rp 


1.43 ±0.03 


1.426 ±0.057 


1.390 ±0.033 


Rjup 


Planet density 


Pp 


0.99 ±0.15 


1.01 ±0.16 


0.93 ±0.12 


gem -3 


Equilibrium Temperature 


T 


1657 ± 55 


1735 ± 34 


1715 ±36 


K 



Notes. (a) As derived from the HARPS/Standard RV measurements. 



change the above results. However, this circumstance was dis- 
cussed in detail in Sects. 3.3.1 and 4.2 o f lLanza etal] (120111) and 
seems to be unlikely. 

In conclusion, our results suggest a smooth, non-violent in- 
ward migration of CoRoT-1 lb during the system formation, in- 
volving a ngular momentu m exchanges with the protoplanetary 
disc (e.g.. lLin et al.ll 1996h . which would not have perturbed the 
planet's primordial alignment. 
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